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We develop a straightforward and quantitative two-step method for spectroscopically classifying galaxies 
from the low signal-to-noise (S/N) optical spectra typical of galaxy redshift surveys. First, using x^-fitting 
of characteristic templates to the object spectrum, we determine the relative contributions of the old stellar 
component, the young stellar component, and various emission line spectra. Then, we classify the galaxy 
by comparing the relative strengths of the components with those of galaxies of known morphological type. 
In particular, we use the ratios of (1) the emission line to absorption line contribution, (2) the young to 
old stellar contribution, and (3) the oxygen to hydrogen emission line contribution. We calibrate and test 
the method using published morphological types for 32 galaxies from the long-slit spectroscopic survey of 
Kennicutt (1992) and for 304 galaxies from a fiber spectroscopic survey of nearby galaxy clusters. From 
an analysis of a sample of long-slit spectra of spiral galaxies in two galaxy clusters, we conclude that the 
majority of the galaxies observed in the fiber survey are sufficiently distant that their spectral classification 
is unaffected by aperture bias. Our spectral classification is consistent with the morphological classification 
to within one type (e.g. E to SO or Sa to Sb) for > 80% of the galaxies. Disagreements between the spectral 
and morphological classifications of the remaining galaxies reflect a divergence in the correspondence between 
spectral and morphological types, rather than a problem with the data or method. 


1. introduction 


in 

,lE}lassification is a basic step toward understanding na- 
tfi^ The classification of galaxies provides a framework 
within which to understand how galaxies form and evolve, 
fundamentals of the morphological classification scheme 
galaxies as proposed by Hubble (1926, 1936) are still in 
pj^e, despite some modifications and additions to the basic 
sequence (cf. de Vaucouleurs 1959; Sandage 1961). The in¬ 
trinsic value of the morphological sequence arises from the 
strong correlations between morphological type and many 
physical characteristics of galaxies (cf. Roberts & Haynes 
1994). Spectral classification, pioneered by Morgan and 
Mayall (1957), provides an alternative to classification based 
on images. The principal advantage of this approach is that 
it is based directly on the physics of the stellar populations 
and the interstellar medium. 

The large volume of data and the need for uniform, 
reproducible, and physically-oriented galaxy classification 
motivate the search for an efficient spectral classification 
method. Such a method, based on the stellar populations 
and emission-line properties of galaxies, must be success¬ 
ful even when applied to the low S/N spectral data typical 
of current galaxy redshift surveys. Spectral classifications 
correlate well with morphological classifications (Morgan & 
Mayall 1957), but also complement them by providing in¬ 
formation not available solely from the morphology. 

We revisit spectral classification armed with digital data, 
computers, and potentially thousands of galaxy spectra. 


Parallel efforts are underway {e.g., Heyl 1994; Connolly et 
al. 1995; Bershady 1995). The procedure described by Heyl 
is potentially quite similar to ours, although his classifica¬ 
tion is based currently only on the galaxies in the Kennicutt 
(1992) sample, which, because of its small size, may not de¬ 
scribe the full range of galaxy properties within any given 
galaxy type. Additionally, our technique has some technical 
advantages, such as the independent determination of emis¬ 
sion and absorption line strengths. The technique developed 
by Connolly et al. primarily utilizes the shape of the contin¬ 
uum (and its sensitivity to recent star formation) to classify 
galaxies. As is evident from their Figure 3, the major dif¬ 
ference between the two principal “eigen-spectra” used to 
classify their galaxies is the continuum behavior blueward 
of 4000 A. This difference arises from the combination of 
the 4000 A break and UV continuum emission from young 
stars. Their technique is optimized for flux-calibrated spec¬ 
tra that extend over a large wavelength region and is best 
suited to spectra that extend blueward of 4000 A in the 
rest frame, such as those of high-redshift galaxies, but gen¬ 
erally applicable to any set of galaxy spectra. Bershady’s 
method is based on galaxy colors, which are observationally 
more efficient for samples of galaxies that are concentrated 
on the sky (such as high redshift or cluster galaxies), and 
bypasses the problem of aperture bias discussed below. On 
the other hand, some degeneracy exists in the colors between 
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the line and continuum contributions, the colors need to be 
K-corrected, and photometric observations are needed in a 
variety of colors. 

In this paper, we describe our method for spectroscopically 
classifying nearby galaxies and examine the correspondence 
between spectral and morphological classification. The prin¬ 
cipal difference between previous work and this study is our 
emphasis on the spectral lines rather than on the contin¬ 
uum shape (although we use some “continuum features”, 
such as the 4000 A break which is actually dominated by 
Ca II H & K absorption). Our philosophy is to keep the 
approach as simple as possible (adopting the fewest fitting 
parameters and the smallest number of classification criteria 
necessary to account for the dominant spectral signatures of 
different morphological types), so that the method is stable 
and can be applied to low S/N data. We fit characteristic 
template spectra to the object spectrum in order to deter¬ 
mine the contributions of different spectral components. We 
then classify each galaxy using three ratios formed from the 
relative contributions of an early-type star (A), a late-type 
star (K), and an emission-line component. Our method can 
be applied uniformly and reproducibly to large numbers of 
the low S/N spectra typical of redshift surveys, which cre¬ 
ates an opportunity to examine the properties of galaxies 
as a function of environment on an unprecedented scale. 
One potential limitation of the method, which we refer to 
as aperture bias, is that the small apertures characteristic 
of fiber or multislit observations may incompletely sample 
the galaxy. Therefore, we investigate the effect of aperture 
bias on our spectral classifications. We conclude that the 
majority are unaffected by aperture bias. 

In §2, we describe the three spectroscopic samples used to 
calibrate and test the method. In §3, we discuss the clas¬ 
sification algorithm, the application to the spectra, and the 
correspondence between spectral and morphological classi¬ 
fication. In §4, we summarize our conclusions. Two appen¬ 
dices contain detailed discussions of the S/N limitations and 
the effect of aperture bias. 

2. THE DATA 

In this section, we describe the three different samples of 
galaxies used to calibrate and test our method of spectral 
classification. The first is a compilation of nearby and well- 
studied galaxies from a spectroscopic survey by Kennicutt 
(1992). We refer to this sample as the Nearby Galaxy (NG) 
sample. These are high S/N spectra obtained by allowing 
the galaxy to drift across the spectrograph slit. This ob¬ 
servational technique produces spectra that are luminosity- 
weighted averages and representative of the entire galaxy. 
Therefore, although these galaxies are typically very nearby, 
the NG sample does not suffer from aperture bias. The sam¬ 
ple’s shortcoming is that its small size prevents us from fully 
sampling the variations in the correspondence between spec¬ 
tral and morphological type for any one type of galaxy. 

The 55 original spectra in the NG sample, kindly provided 
by R. Kennicutt, were obtained mostly at the Steward 2.3m 
telescope and cover the wavelength region from 3650A to 
7100A. The data are flux calibrated, but are not corrected 
for either Galactic or internal reddening (for details see Ken¬ 
nicutt 1992). The spectra include all of the major optical 


lines of galaxy spectra, except the Ca triplet at about 8500A 
(cf. Figure 1 for example spectra). The spectral resolution 
is between 5 and 8 A, which is comparable to that typically 
obtained in redshift surveys. 

For the purpose of calibrating and testing our classifica¬ 
tion algorithm against the NG sample, we study only the 32 
“normal” galaxies, where “normal” denotes a lack of broad 
emission lines. These are the first 27 galaxies in the Kenni¬ 
cutt catalog (1992; see his Table 1), plus NGC 2798, NGC 
3034, NGC 3077, NGC 3310, and NGC 5195. Morphological 
types for these galaxies are drawn from Kennicutt’s Table 
1 and are from either the Revised Shapley-Ames Catalog 
(Sandage & Tammann 1981), the Second Reference Cata¬ 
log of Bright Galaxies (de Vaucouleurs, de Vaucouleurs, & 
Corwin 1976), Huchra (1977), or Kennicutt (1992). Our re¬ 
binning of these morphological types into a coarser grid is 
outlined in Table 1. Based on the published morphological 
classifications and our rebinning of types, the NG sample 
contains 4 E’s, 3 SO’s, 5 Sa’s, 5 Sb’s, 9 Sc’s, and 6 Irr’s. 


Table 1 : Our Rebinning of Galaxy Types 


Our Type 

Standard Type 

T-Type 

E 

E(dwarfs and giants) 

-5,-4 

SO 

L, SO 

-3,-2,-1 

Sa 

S0/a,Sa 

0,1 

Sb 

Sab,Sb 

2,3 

Sc 

Sbc,Sc 

4,5 

Irr 

Scd,Sd,Sm,Irr 

6-10 


Our second and much larger sample consists of 304 new 
fiber spectra of galaxies with published morphological clas¬ 
sifications from Dressier (1980) and Richter, Materne, & 
Huchtmeier (1987) in the fields of nearby galaxy clusters. 
We refer to this sample as the Cluster (CL) sample. This 
sample is a factor of 10 larger than the NG sample and 
thus allows us to statistically examine the correspondence 
between spectral and morphological classification for galax¬ 
ies of each type. This sample’s potential shortcomings are 
that classification errors could result from either the lim¬ 
ited fiber aperture size (aperture bias) or the low S/N of 
the spectra and that the spectral coverage (3300A to 6200 
A) does not extend redward to Ha. However, because these 
data are very similar in quality to those of redshift surveys 
of large-scale structure, this sample provides a fair test of 
our spectral classification algorithm. 

The CL sample spectra, which are part of a larger survey 
of the internal dynamics of rich clusters of galaxies (Zablud- 
off 1995), were obtained at the Las Campanas Observatory 
(LCO) with the 2.5m telescope and multifiber spectrograph 
(Shectman et al. 1992). The recessional velocities of galax¬ 
ies in the survey are mostly between 3000 and 50000 km 
s^^. The fiber aperture of 3 arcsec corresponds to a range 
of physical sizes between ~ 0.6 and 10 kpc, respectively (Hq 
= 75 km s“^ Mpc“^). The galaxies have B-band apparent 
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Figure 1: Characteristic spectra from the Kennicutt (1992) sample of normal galaxies of types E to Irr. The wavelength 
coverage of each panel is from 3600 to 7100A. Even among adjacent types there are noticeable spectral differences. 


magnitudes between ~ 15 and 18.5, and the total exposure 
times are between 2 and 6 hours. These data are not flux 
calibrated, which is typical of the redshift survey data to 
which our classification method might be applied. The lack 
of flux calibration prevents us from using the continuum 
shape as a classification criterion. 

The third sample is a random subsample of 18 new long- 
slit spectra of late-type galaxies from a large sample ob¬ 
served as part of a Tully-Fisher survey of galaxies in clusters 
with recessional velocities of cz ~ 10000 km s^^ (Willick 
1995). We refer to this sample as the aperture bias (AB) 
sample. These spectra were obtained with the Palomar 5-m 
telescope and double spectrograph. We use these data to 
quantify the effect of aperture bias on classifications made 
with small aperture spectra. The spectral coverage is from 
3850 to 4700 A and from 6450 to 7050 A, with a resolution 


between 1 and 2 A. The 1 arcsec wide slit was oriented along 
the major axis of the target galaxy. The exposure times are 
typically between 10 and 30 min. We have accompanying 
CCD images obtained at the Palomar 1.5m telescope from 
which we determine the morphological classifications and 
major axis orientation. These spectra are typically of lower 
S/N than the NG sample, but of higher S/N than the CL 
sample. To test how particular spectral features vary with 
galactocentric radius, we use the blue spectra to examine 
absorption line (Ca II H & K, and G-band) characteristics 
and the red to examine emission line (Ha) characteristics. 
These data are also not flux calibrated. An examination of 
aperture bias using this sample is presented in Appendix B. 

The spectra in the CL and AB samples were reduced with 
IRAF in a standard manner for multifiber and long-slit data, 
respectively. The reductions consisted of bias subtraction. 
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flat flelding, wavelength calibration, extraction, and sky sub¬ 
traction (the last two steps are reversed for the long-slit 
data)^. 

3. CLASSIFYING GALAXIES 

Our approach to spectral classification assumes only that 
galaxy spectra are a combination of (1) an old stellar spec¬ 
trum, which we characterize as a K-star spectrum, (2) a 
young stellar spectrum, which we characterize as an A-star 
spectrum, and (3) an emission line spectrum. In the simplest 
terms, our classification procedure consists of measuring the 
relative contributions of these components to the observed 
galaxy spectrum and then comparing their relative strengths 
to those in the spectra of galaxies of known morphological 
type. We base this approach on the most striking spectral 
differences among the galaxies in the NG sample. 

We use x^-fltting to determine which linear combination 
of template spectra best reproduces the observed spectrum. 
Our set of templates consists of a composite^ spectrum of a 
K star (characteristic of a stellar population > 5 Gyr old), a 
composite spectrum of an A star (characteristic of a stellar 
population < few Gyr old), and an H II region (characteristic 
of ongoing star formation) divided into four separate tem¬ 
plate spectra: (1) the [O II]AA3726,3729 lines, (2) the II/3 
line, (3) the [O III] AA4959,5007 lines, and (4) the Ha line. 
We excise from each composite stellar template the principal 
characteristic lines from the other template to construct two 
nearly “orthogonal” template spectra (the correlation be¬ 
tween the final two stellar templates is 0.13). The modified 
K-star composite template includes all of the stellar absorp¬ 
tion lines except for the hydrogen Balmer lines, while the 
composite A-star template includes only the Balmer lines. 
The H II region template is taken from an Irr galaxy (NGC 
1569) in the Kennicutt sample. The quality of the GL data 
does not warrant the use of more detailed template spec¬ 
tra {e.g., separate H & K and Mg 2 templates). We justify 
the choice of these templates over a mathematically orthog¬ 
onal set derived from the data themselves (such as done by 
Connolly et al. ) by noting that for our low S/N spectra 
much of the variance is due to the noise, typically about one- 
half of the variance is accounted for by our best fit, so that 
a principal component analysis of the CL sample will not 
efficiently translate the mathematical properties of the sam¬ 
ple into the physical properties of the underlying parameter 
space. A principal component analysis, such as that of Con¬ 
nolly et al. , on a large sample of high-S/N spectra would 
test the assumptions involved in our template selection. 

We prepare each object and composite template spectrum 
before applying the x^-Atting. First, we continuum subtract 
and standardize the template spectra to have a mean » 0 
“counts” (for technical reasons) and a small maximum 

^IRAF is distributed by the National Optical Astronomy Observato¬ 
ries, which are operated by the Association of Universities for Research 
in Astronomy, Inc. (AURA) under cooperative agreement with the Na¬ 
tional Science Foundation. 

^The composite absorption-line templates are made from serendipi¬ 
tous observations of K-like and A-like stars with the Las Campanas 
multifiber system (S. Shectman 1995, private communication). 


deviation from this mean (fixed, but arbitrarily chosen). 
The particular numbers in this procedure have no specific 
meaning, other than to insure the uniformity of the data. 
Second, the galaxy spectra are continuum subtracted and 
standardized to have an rms dispersion equal to the maxi¬ 
mum deviation from the mean in the template spectra. Gon- 
tinuum subtraction, rather than division, is used because (1) 
continuum variations between objects over a limited wave¬ 
length region are primarily due to intrinsic brightness and 
color changes rather than detector sensitivity and (2) the 
low continuum levels in some objects would introduce sig¬ 
nificant noise. The spectra are rescaled to a fixed rms so 
that there is some standardization among template ampli¬ 
tude fits, but this is not a necessary step. The spectra are 
now ready for the x^-fitter. 

We At our templates with the x^-minimization fitting al¬ 
gorithm kindly provided by H.-W. Rix, which was developed 
by Rix & White (1992) to study detailed kinematics in ellip¬ 
tical and SO galaxies. We At the absorption line components 
flrst, by random walking through a parameter space that in¬ 
cludes the recessional velocity, the Gaussian line-broadening 
function, and the relative contributions of the two absorp¬ 
tion line templates. The best At is then subtracted from 
the original spectrum. We then At the emission lines to the 
residual spectrum, allowing the broadening function to dif¬ 
fer from that used for the absorption lines. (The emission 
line broadening, which may arise from a few hot regions of 
narrow (~ 10 km s^^) line width, may difler from the stel¬ 
lar velocity dispersion (the terminology here includes bulk 
rotation) characterized by the absorption line broadening 
function.) The best-flt spectrum of a galaxy is then the 
sum of the absorption-line and emission-line flts created by 
minimizing x^ on the two separate passes. Because of the 
two step fltting procedure, there is no possibility that emis¬ 
sion and absorption lines would be constructed to cancel 
each other. 

The x^'fitting routine quantifles the contribution from 
each template, which we call template amplitudes, to the 
best-flt model spectrum for each galaxy. The template am¬ 
plitudes are > 0. We calculate a goodness-of-flt criterion by 
dividing the rms deviation in the original spectrum by the 
rms deviation in the residual spectrum. The At is acceptable 
if the rms residual is at most 70% of the original (a value 
determined from visual inspection to return reliable, S/N > 
3, spectral features). Internal errors can be estimated from 
the range of models with a value of x^ that is acceptable 
at some confldence level. One could propagate this uncer¬ 
tainty through the classiflcation procedure to determine an 
internal estimate of the classiflcation uncertainty. However, 
this is a complex procedure that relies on several estimated 
uncertainties and so we prefer to rely upon the external un¬ 
certainty estimates that are discussed further below. Figure 
2 shows the At for a random Sb galaxy in the Kennicutt 
sample (NGG 3147), whose rms residual is 65% of the rms 
in the original spectrum. 

Several important aspects of the fltting routine are appar¬ 
ent in Figure 2. First, the continuum subtraction removes 
all structure on scales > 200 A and in doing so does not 
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Figure 2: The original, continuum subtracted, and model spectra of a representative Sb galaxy (NGC 3147) from the NG 
sample. In the right panels, the dashed lines trace the fitted spectra. 


alter any spectral line feature. Second, the model spectrum 
has almost no structure for A > 5500 A because our absorp¬ 
tion line templates do not extend beyond 5500 A. Third, 
the strengths of the absorption lines are somewhat under¬ 
estimated. The lower right panel shows the fit to the Ca 
II H and K lines in greater detail. The slight discrepancy 
between the data and the fit is in part the result of the in¬ 
version of the ratio of H to K strengths between the K-star 
template and the object and of the imperfect matching of 
the width of the features. The latter may be the result of 
a non-Gaussian line profile, as one might expect from disk 
rotation in an Sb galaxy. Lastly, we do not include an [N II] 
template, so the fit around Ha is often distorted in an at¬ 
tempt to account for some of the [N II] flux with the Ha line. 
Despite these small flaws in the fit, we classify the galaxy 
as an Sb (as described in detail below), which matches its 
morphological classification (Kennicutt 1992). 

Our spectral classification is based on three ratios: (1) the 
sum of the emission-line template amplitudes to the sum of 
the absorption-line template amplitudes, (2) the A-star tem¬ 
plate amplitude to the K-star template amplitude, and (3) 
the sum of the oxygen emission-line template amplitudes to 
the sum of the hydrogen emission-line template amplitudes. 


These ratios were chosen because they isolate the most no¬ 
ticeable differences among the galaxy spectra in Figure 1 
and represent physical components (ongoing star-formation, 
recent star formation, the dominance of an older stellar pop¬ 
ulation, and the excitation of the gas) in the spectra. 

The logarithms of the three ratios form a 3-dimensional 
space in which the galaxies are distributed. Our basic as¬ 
sumption is that galaxies of each morphological class in¬ 
habit a reasonably distinct region of this space — an as¬ 
sumption that is supported by the range of spectral char¬ 
acteristics seen in Figure 1 and by the results of Morgan & 
Mayall (1957). We define that region for a particular class of 
galaxy by calculating the positional centroid and the disper¬ 
sion along each of the three coordinate axes for galaxies of 
that morphological type. Regions are defined for E, SO, Sa, 
Sb, Sc, and Irr galaxies using the published morphological 
classifications. A galaxy is spectrally classified by assigning 
it to the nearest group in this ratio space. The distance 
to any particular group is calculated by taking the distance 
along each of the three coordinate axis between the galaxy 
and the centroid of the group, dividing those distances by 
the dispersion of the group along that axis, and then adding 
the three components in quadrature. The distance along 
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a particular axis is not used if that ratio is ill-defined; for 
example, if a galaxy does not have emission lines, then the 
ratio of oxygen to hydrogen emission is not used and that 
galaxy is classified within the 2-D space defined by the other 
two ratios. A projection of these regions is provided in §3.2 
with the discussion of the classification of the CL sample. 

Several technical details require further clarification. 
First, we use the NG and CL samples to self-calibrate the 
classification regions. Second, we use template amplitudes 
rather than equivalent widths because the continuum, which 
is often vary weak, would introduce large uncertainties into 
equivalent widths. Third, we use the logarithm of the tem¬ 
plate amplitude ratios to linearize differences among ratios. 
Fourth, when a template amplitude is less than 0.001, the 
template contribution is set at 0.001. This limit avoids di¬ 
vergent logarithms. Fifth, we determine whether a spectrum 
has emission lines by checking whether the largest emission 
line template amplitude is at least 0.20 (which translates 
typically into a S/N > 3 in the peak of the line) or whether 
the median emission line template amplitude is greater than 
0.001 (which implies that the program found two or more 
candidate emission lines corresponding to the absorption 
line redshift). Finally, we examined a variety of different 
combinations of template amplitudes ratios before deciding 
on the choices outlined above. 

a) the nearby galaxy (ng) sample 

The NG sample, the simplest of the samples to analyze, 
provides the cleanest test of the method. After applying 
our spectral classification algorithm, 26 out of the 32 galax¬ 
ies are classified to within one morphological type (81%). 
“Misclassifications” can arise from (1) erroneous morpho¬ 
logical classifications, (2) errors in the fitting procedure, (3) 
errors in the classification procedure, and (4) scatter in the 
relationship between the morphological and spectral prop¬ 
erties of galaxies. We neglect the first because it is beyond 
our control. The other sources of “misclassification” are the 
focus of the remainder of this paper. 

There are several possible sources of errors in the fitting 
procedure. First, the templates have a slightly different 
spectral resolution than the NG data. Second, for simplicity 
we have not included a [N II] template, so the Ha fits are 
sometimes distorted in an attempt to account for some of 
the [N II] emission. Third, the velocity broadening function 
is taken to be a Gaussian, which may lead to errors in the 
fit of galaxies with significant rotation. Fourth, we make no 
allowance for differences between the velocity widths of per¬ 
mitted and forbidden emission lines. Finally, we assume im¬ 
plicitly that all the absorption lines characteristic of the old 
population {e.g., Ca II H & K, Mg 2 , Fe) are characterized by 
a single template spectrum. For all of these reasons, there 
is likely to be some template mismatch. We have chosen to 
minimize the number of templates and fitting functions in 
order to retain the stability of the procedure, rather than 
attempt to address some of these possible sources of uncer¬ 
tainty. These issues are mostly irrelevant for the low-S/N 
data of the CL sample because other sources of uncertainty 
dominate (cf. §3.2 and Appendices). 

Uncertainties in the classification procedure are intro¬ 
duced by the small size of the NG sample. There are typi¬ 


cally only three to five galaxies in each type bin. Therefore, 
the delineation of the regions of ratio space occupied by each 
morphological type is highly uncertain, and one peculiar ob¬ 
ject easily skews the classifications. Finally, we note that the 
same data are used both to classify and test the procedure. 
This approach may artificially elevate the success rate in 
such a small sample. These problems are less severe for the 
CL sample because it is a factor of 10 larger. 

After inspection of individual fits and classifications for 
the NG sample, we conclude that the disagreements between 
the spectral and morphological classifications do not gener¬ 
ally reflect problems in the quality of the data or the algo¬ 
rithm, but instead reveal a break in the detailed correspon¬ 
dence between spectral and morphological properties. Some 
of these galaxies are known to be morphologically peculiar. 
These galaxies are included in the NG sample because we 
chose to select our “normal” sample based solely on spec¬ 
troscopic criteria, which we could apply to any sample of 
galaxies. The NG sample includes three objects that Ken- 
nicutt classifies as peculiar, one object that RSA classifies 
as amorphous, and the companion to M 51. Two of our six 
“misclassifications” are from this peculiar subsample. 

The remainder of the “misclassifications” in the NG sam¬ 
ple arise from scatter in the spectral characteristics within 
each morphological type. The regions occupied by galaxies 
of different morphologies in the 3-D ratio space overlap. It is 
therefore unavoidable that galaxies from one type will spill 
over into the region occupied by another type. We conclude, 
based on our inspection of the observed and model spectra, 
that the dispersion in spectral properties is real and not an 
artifact of the fitting procedure. 

b) the cluster (cl) sample 

We now calibrate and test the classification method with 
the larger CL sample. Because the CL spectra do not extend 
beyond 6200 A, we do not use the Ha emission line template. 
Otherwise, the templates are the same as those used for the 
NG sample. 

We find that the spectral and morphological types are con¬ 
sistent for 253 of the 304 sample galaxies (83%). This frac¬ 
tion is almost identical to that found for the much smaller 
NG sample. The centroids and standard deviations of the 
three template amplitude ratios (as derived from the NG 
and CL samples) are listed in Table 2 for each morpholog¬ 
ical type, and a projection of the regions occupied by the 
galaxies is shown in Figure 3. The ratios differ between the 
NG and CL samples because of the inclusion of Ha in the 
analysis of the NG sample. 

The oxygen to hydrogen emission line ratio has no dis¬ 
criminatory power for this sample (the same classification 
“success” rate is achieved without this ratio). Nevertheless, 
we retain it in the classification algorithm because it does 
vary among the latest types in the NG sample (cf. Figure 1) 
and would probably increase the number of successful clas¬ 
sifications if the CL sample had a larger fraction of late-type 
spirals. For comparison, using only the emission to absorp¬ 
tion template amplitude ratio produces a “misclassification” 
fraction of 0.25 for the CL sample and using only the A-star 
to K-star ratio results in a “misclassification” fraction of 
0.61. 
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Table 2 : Spectral Classification Ratios from NG and CL Sample Galaxies 


Type 

log(Em/Ab) 

log(K Star/A Star) 

log(0 Em/H Em) 

Number of Galaxies 

NG Sample 
(without Ha) 






E 

-2.52 

± 0.10 

-2.08 ± 1.07 

— 

4 

SO 

-2.52 

± 0.10 

-0.75 ± 0.28 

— 

3 

Sa 

-1.85 

± 1.51 

-0.92 ± 0.65 

-0.07 ± 0.10 

5 

Sb 

-0.49 

± 1.21 

-0.69 ± 0.88 

0.75 ± 0.92 

5 

Sc 

1.13 

± 0.69 

-0.59 ± 0.68 

0.16 ± 0.23 

9 

Irr 

0.82 

± 1.79 

-0.44 ± 0.78 

0.42 ± 0.22 

6 

CL Sample 






E 

-2.42 

± 0.61 

-1.52 ± 1.10 

0.24 ± 0.23 

52 

SO 

-2.38 

± 0.82 

-1.11 ± 0.85 

0.55 ± 1.68 

151 

Sa 

-1.83 

± 1.24 

-1.06 ± 0.98 

0.32 ± 1.02 

40 

Sb 

-1.20 

± 1.53 

-0.92 ± 1.05 

1.32 ± 1.65 

30 

Sc 

-0.23 

± 1.06 

-0.50 ± 0.87 

0.31 ± 0.92 

26 

Irr 

0.11 

± 0.19 

-0.48 ± 1.31 

0.36 ± 0.27 

5 


In Figure 4, we plot the fraction of galaxies whose spectral 
and morphological classifications agree to within one type 
as a function of morphological type. The Figure illustrates 
that the CL sample is heavily biased toward elliptical and 
lenticular galaxies. The agreement between the spectral and 
morphological types is significantly greater at either end of 
the type range than in the middle because early-type spi¬ 
rals exhibit the largest variations in spectral properties (as 
evident in Figure 3), especially when Ha, which is brighter 
and less obscured by dust than H/3, is excluded from the 
analysis. These intermediate-type galaxies have also been 
the most difficult to spectrally classify in previous studies 
(Connolly et al. and Bershady) and are the most difficult to 
classify visually. The similarity in the success rates between 
the NG and CL samples suggests that neither low S/N nor 
aperture bias has seriously affected our classification of the 
CL galaxies. 

We now test the extent of the effects of low S/N and aper¬ 
ture bias on the classification of the CL sample. In the 
left panel of Figure 5, we plot a histogram of the classifi¬ 
cation success as a function of the goodness-of-fit criterion 
discussed previously (larger numbers indicate a better fit). 
This criterion reflects both the quality of the original spec¬ 
trum (the S/N) and the success of the fitting algorithm, 
which in turn depends on such factors as template and ve¬ 
locity broadening mismatch. There is no evidence within the 
range spanned by the data that there is any dependence of 
the classification algorithm on this ratio. This result demon¬ 
strates that many “misclassiflcations” have not been intro¬ 
duced by the inclusion of poor quality spectra. The effect 


of noise on the fitting procedure is mostly to introduce spu¬ 
rious emission-line detections, which we minimized by the 
judicious selection of our emission-line criteria discussed at 
the beginning of §3. In Appendix A, we discuss adding noise 
to the NG spectra and demonstrate that the classification 
algorithm is reliable for S/N >3. 

Because the galaxies in the CL sample are all observed 
with the same size aperture, the spectrum of a higher red- 
shift galaxy includes light from a larger physical radius than 
does the spectrum of a lower redshift galaxy. To examine the 
effect of aperture bias on the CL sample, we plot a histogram 
of the classification success as a function of redshift in the 
right panel of Figure 5. The success rate is fairly constant 
across the range of redshifts and the misclassiflcations do 
not appear to depend on morphological type. Aperture bias 
does not systematically affect the classification of galaxies in 
the CL sample, and there appears to be no justification for 
excluding galaxies with 2 ; < 0.1 from this type of classifica¬ 
tion scheme. However, this statistical result does not prove 
that aperture bias is unimportant in the classification of any 
individual galaxy; in particular, we And that aperture bias 
can significantly distort some galaxy spectra (cf. Appendix 
B). 

In Figure 6, we plot a comparison between the morpho¬ 
logical and spectral classifications for the 51 galaxies that 
were “misclassifled”. Note that there are galaxies that are 
widely discrepant; for example, two morphological E’s are 
spectrally classified as Irr’s (see Figure 7 for the spectra of 
the three galaxies with the most discrepant spectral and 
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log(Emission/Absorption Template Amplitudes) 

Figure 3: The spectral classification criteria derived from the 
CL sample projected onto the plane of emission/absorption 
and early/late template amplitudes. The ellipses are cen¬ 
tered at the mean position for galaxies of each morphological 
type and are scaled to match the standard deviation along 
each axis of the galaxies within that type. 



S/N Velocity (km/sec) 


Figure 5: The classification success rates for the CL sample 
as a function of goodness-of-fit (as described in text) and 
redshift. The shaded regions represent the successful clas¬ 
sifications. The rightmost bin in each panel includes all of 
the galaxies with characteristics that lie to the right of this 
bin. The fraction of successful classifications is given above 
each bin. In the upper portion of the right panel, we plot 
the morphological types of the galaxies misclassified in the 
corresponding bin. 



Figure 4: The fraction for which the morphological and spec¬ 
tral classifications agree to within one type as a function of 
morphological type for the CL sample. The shaded areas 
represent the successful fraction (the numerical fractions are 
also given). 
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Figure 6: A comparison between the spectral versus morpho¬ 
logical classification for the 51 galaxies from the CL sam¬ 
ple for which the spectral type varied by more than one 
type from the morphological classification. The filled cir¬ 
cles represent the morphological type from Dressier (1980) 
or Richter, Materne, & Hucthmeier (1987), and the open 
circles represent our spectral classification. 
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Figure 7: The spectra of the three galaxies in the CL sample 
with the most discrepant spectral and morphological types. 
Major spectral lines are labeled. The morphological type is 
shown to the left of the arrow and the spectral type to the 
right. The spectra have been smoothed by a factor of 4. 

morphological types). In future applications, it might be 
possible prevent the classification of a typical E as an Irr 
by introducing a magnitude criterion (or the velocity broad¬ 
ening determined with the x^-fitter). There are a total of 
13 morphological E’s or SO’s that are classified as spec¬ 
tral type Sb or later. The algorithm has “found” emission 
lines in all 13, although only eight have obvious emission 
lines (the remainder are marginal line detections). If the 
five galaxies without obvious emission are all spurious de¬ 
tections, the emission-line false detection rate is 2% (5 out of 
203) among the E’s and SO’s. With our data, we cannot de¬ 
termine whether the eight galaxies with obvious emission are 
morphologically misclassified as early-types or are in fact el¬ 
lipticals with some ongoing star formation (none have broad 
emission lines). The remaining widely discrepant galaxies 
are 12 late type galaxies (Sb or later) with strong H and 
K lines and no emission lines that we spectrally classify as 
either E’s or SO’s. It is unclear whether the spectral classifi¬ 
cation of these spirals is affected by aperture bias or reveals 
old disks with very low star formation rates. 

These results, and those for the NG sample, indicate that 
the spectral and morphological classifications disagree by 
more than one type for only about one-fifth of the galaxies. 
As mentioned in §3.1, this inconsistency can arise for a va¬ 
riety of reasons. The interesting scientific issue is to what 
degree these mismatches result from a break in the corre¬ 
lation between morphological and spectral characteristics. 
Resolving this issue requires a detailed knowledge of the 
contribution of the possible sources of error. Based on the 
arguments presented in the Appendices and the very sim¬ 


ilar misclassification fractions in the NG and CL samples, 
we conclude that most of the misclassifications arise from 
real scatter in the spectral properties of galaxies within each 
morphological bin. Some of that scatter is relatively small as 
evidenced by the overlap seen in Figure 3 (only 24% of the 
CL misclassifications are due to mismatches of more than 
2 types). However, there is a small fraction of galaxies for 
which the spectral and morphological types are grossly dif¬ 
ferent, such as the morphological E’s with emission lines. 
Such objects require further study and demonstrate why 
spectral and morphological classifications are complemen¬ 
tary. 

Lastly, we present one more test of the classification al¬ 
gorithm. We use the classification calibration derived from 
the CL sample to classify the NG sample. We recalculate 
the various template amplitude ratios for the NG sample 
excluding Ha, since Ha is not used in the CL sample. Ap¬ 
plying the CL criteria we correctly classify 91% of the NG 
sample. This success rate is higher than that attained with 
self-calibration (cf. §3.1). This result partially addresses the 
concerns that the lack of flux calibration severely distorted 
the template amplitude ratios and that the success rate has 
been elevated by using self-calibration. We also present in 
Table 2 the calibration derived from the NG sample exclud¬ 
ing Ha to enable a direct comparison of the derived criteria 
from the two samples. The values from the NG sample have 
larger external uncertainties due to the small number of ob¬ 
jects per bin. We conclude by stressing that the applicability 
of the classification criteria derived for the CL to other sam¬ 
ples has not been fully investigated. Any application of this 
technique to other samples must be self-tested. 

4. SUMMARY 

We develop and test a method of spectroscopically classi¬ 
fying nearby galaxies. The strengths of the method are that 
(1) for each galaxy, we obtain an estimate of the relative 
contributions of the populations of young and old stars, (2) 
unlike morphological classification, the method is indepen¬ 
dent of spatial resolution, (3) the method can be applied 
uniformly to existing data for thousands of galaxies, (4) the 
classification uncertainties are quantifiable through Monte- 
Carlo simulations and comparison with morphological clas¬ 
sifications, and (5) the method is independent of continuum 
shape and so does not require fiuxed spectral data. There 
is a strong correspondence between spectral classification 
and visual morphological classification (> 80% of galaxies 
can be spectrally typed to within one morphological type). 
The spectral types might correspond even more strongly to 
other galaxy properties that are known to correlate with 
Hubble type {e.g., HI content, environment, or luminosity), 
and an examination of those correlations can provide greater 
insight into galaxy evolution. In future applications, the in¬ 
dependence of our method on continuum shape might also 
be utilized to test the method against the other classification 
methods that depend on continuum shape. 

The principal applications of spectral classification algo¬ 
rithms are to high-redshift surveys, where the limited spatial 
resolution prohibits visual classification, and to large red- 
shift surveys. Whether one uses colors, continuum shape, or 
spectral line properties as the basis of spectral classification 
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depends primarily on the type of data that is available. It is 
now possible to determine the local density-type, luminosity- 
type, and correlation function-type relations from large red- 
shift surveys without visually classifying tens or even hun¬ 
dreds of thousands of galaxies. In addition, these same re¬ 
lationships can be examined in terms of spectral character¬ 
istics, such as the presence of a young stellar population. 

We define and test our algorithm with three samples. The 
first is a small, but otherwise ideal sample (NG) that does 
not suffer from aperture bias and has high S/N. The re¬ 
sults are used to illustrate the success of the basic classifi¬ 
cation procedure. The second sample (CL) consists of low 
S/N fiber spectra of galaxies of known morphological type 
taken from a redshift survey with a multifiber spectrograph. 
These spectra are of similar quality to those obtained from 
large redshift surveys {e.g., Huchra et al. 1992; Shectman 
et al. 1992). We find that despite an important limitation 
of the data, the lack of spectral coverage at the Ha line, our 
spectral classification matches the morphological type for > 
80% of the CL galaxies (the same as the fraction “correctly” 
classified in the NG sample and consistent with the results 
of Morgan & Mayall 1957). With these two samples and a 
third sample (AB) of long-slit spectra, we quantify the ef¬ 
fects of low S/N and of aperture bias on the classifications. 
We find that the algorithm performs well in the low S/N 
regime (S/N > 3) and that aperture bias does not seriously 
affect the fiber data. An amalgamation of continuum and 
line spectral classification methods should strengthen the 
classification procedure (although fluxed spectra are usually 
not available from large redshift surveys). Finally, we stress 
that although the classification method is general, the am¬ 
plitude ratios and the particular classification criteria we 
derive are not, because the CL data are not flux calibrated. 

The correspondence between spectral and morphological 
classifications places a strong constraint on models of galaxy 
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formation and evolution, because the absence of large popu¬ 
lations of very young E’s or very old S’s is still unexplained. 
The well-known decrease in star formation activity along the 
Hubble sequence is clearly visible in the second and third 
columns of Table 2. In the progression from early- to late- 
type galaxies, both the mean emission to absorption line 
ratio and the A star to K star ratio systematically increase 
(although the scatter is large). 

The comparison between the spectral and morphological 
classifications also reveals some unusual galaxies. For ex¬ 
ample, the samples include ellipticals with line emission and 
late-type spirals in which star formation appears to have 
ceased. In total, about 10% of E’s and possibly a similar 
number of spirals (depending on the degree of aperture bias) 
have more in common with galaxies of the opposite spectral 
type than with their own. Whether these unusual objects 
result from recent galaxy-galaxy interactions or represent 
the respective tails of the evolutionary histories of elliptical 
and spiral galaxies is unknown. This issue may be resolved 
soon, when the use of automated spectral and morpholog¬ 
ical classification (Abraham et al. 1994; Han 1995) makes 
possible a detailed study of the properties and environments 
of thousands of galaxies. 
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Appendix A: Spectral Classification in the Low 
S/N Regime 

Most of the data to which our spectral classification method 
might be applied are low (< 10) S/N spectra obtained orig¬ 
inally for redshift measurements. We now test the limit in 
S/N for which we can obtain reliable classifications by de¬ 
grading the NG spectra to even lower signal-to-noise than is 
characteristic of the CL sample. 

We create a series of noise-enhanced versions of the NG spec¬ 
tra by adding Gaussian uncorrelated noise to each pixel and 
degrading the typical S/N to 3.5, 2.3, and 1.8 in three dif¬ 
ferent tests. These S/N values are based on the ratio of the 
depth of the Ca II absorption lines to the noise — not the 
ratio of the integrated signal in the features to the noise. 
We prefer to use the feature depth rather than its equiva¬ 
lent width because many of the spectra have low continuum 
levels which introduce large uncertainties into the equiva¬ 
lent width. The fitting and classification are then redone by 
using the regions in ratio space defined by the original NG 
data. Decreasing the S/N to 2.3 is sufficient to boost the 
fraction of spectra for which the spectral and morphological 
classifications disagree by more than one type to above 30%, 
which roughly corresponds to a la uncertainty per type in 



Figure 8: The degraded (S/N = 2.3), smoothed, and model 
spectra for NGC 3147. This plot represents the level of noise 
at which significant (> 30%) misclassification occurs. The 
spurious emission lines that produce the misclassification 
can be seen in the model spectrum. 


We quantify how the spectrum changes with noise by cal¬ 
culating the cross-correlation of the original and noise- 
enhanced spectra. The mean correlation value, for the S/N 


= 2.3 case, is 0.59, and so the noise has significantly altered 
the spectrum. Figure 8 shows the noise-enhanced version 
of the spectrum of NGC 3147 shown in Figure 1. We also 
plot a smoothed version to demonstrate that there is still 
significant signal in the spectrum, and then plot the model 
fit. In this case, the problem is that noise spikes are easily 
identified as [O II], H/3, and [O III] emission lines (the Ha 
emission line is not fit because the wavelength range for the 
CL sample does not extend that far redward). We conclude 
that S/N > 3 (peak to noise) is necessary to maintain the 
standard 80% correspondence between morphological and 
spectral properties (cf. Morgan & Mayall 1957; our NG and 
CL results). 

Appendix B: Aperture Bias 


Aperture bias is a potentially serious problem for spec¬ 
tral classification based on multi-object observations. Fiber 
apertures are typically about 2 arcsec (although the LCO 
spectrograph has 3 arcsec apertures), and multislits are usu¬ 
ally also small. At recessional velocities as large as 30000 
km s“^ (z = 0.1), a standard 2 arcsec aperture samples a 
radius of only 2 kpc. 

Our investigation of aperture bias focuses on the AB sam¬ 
ple. This sample includes only spirals, which are generally 
more susceptible to aperture bias than ellipticals, because 
spiral spectral properties are typically less homogeneous as 
a function of radius {e.g., bulge, disk, HII regions, and dust). 
With long-slit data, we can vary the effective aperture and 
thus the area of the disk that we use to classify the galaxy. 

The AB sample consists of 18 galaxies from two galaxy clus¬ 
ters, Abell 2199 (z = 0.03) and Abell 2657 (z = 0.04). For 
each galaxy, there is a blue and red spectrum (cf. §2). With 
these spectra, we quantify the variations in the emission 
and absorption lines as a function of radius. Because of 
the restricted wavelength range of these data, especially the 
absence of the region from 4700A to 5100A, we cannot ap¬ 
ply the classification algorithm directly. Instead, we mea¬ 
sure the change in each spectrum as the effective aperture 
is varied. Each spectrum is reduced in a standard manner 
through the sky subtraction step. We then extract differ¬ 
ent width apertures that trace the galaxy spectrum. The 
aperture width is incremented by 2 pixels for each extrac¬ 
tion until only noise is being added to the total spectrum 
(as defined by the first net decrease in the total counts in 
the spectrum). For each galaxy, we cross-correlate each of 
the extracted spectra with the widest-aperture spectrum to 
estimate the difference between the two. 

Ideally, we would have spectra that fully sample the galaxy, 
rather than long slit spectra. For face-on galaxies, one could 
multiply the ffux observed at radius R within the slit by R 
to get the ffux in an annulus. However, the AB galaxies 
(by virtue of being candidates from a Tully-Fisher study) 
are at large inclinations to the line-of-sight. Any transfor¬ 
mation from the ffux within the slit to total ffux requires 
not only azimuthal corrections, but inclination corrections 
requiring assumptions about the optical depth of the disk. 
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Such corrections to the noisy spectral data at large radius 
are dubious. Therefore, we make no correction to the data 
as a function of radius, but note that the contributions to 
the flux from large radii may be underrepresented. 

Differences among the extracted spectra of a galaxy can 
arise from changes in the S/N. To account for a decrease 
in the cross-correlation value due to noise, we first construct 
a simulated narrow aperture spectrum by adopting the refer¬ 
ence (wide aperture) spectrum as the true spectrum within 
the narrow aperture {i.e., we assume that the signal does 
not change across the entire aperture and that the reference 
spectrum is noiseless). To correct the signal strength for the 
change in aperture size, we divide the simulated spectrum 
by the ratio of the total counts in the reference to those in 
the observed narrow aperture spectrum. We subtract a sim¬ 
ulated sky (with a mean of 0 counts and a noise level given by 
the readnoise Xi/width of the aperture in pixels) and add 
Poisson noise to the signal. The cross-correlation between 
the simulated spectrum (5) and the reference spectrum (R) 
then quantifies the degradation of the cross-correlation due 
to a change in the S/N. The cross-correlation coefficient for 
the observed narrow aperture spectrum (TV) and the ref¬ 
erence is standardized using the previous estimate of the 
degradation due solely to noise. The corrected correlation 
values are 


~ EiSiRiEi^iNi’ 


where the sum over i represents the sum over individual 
pixel values and each spectrum has had a continuum fit sub¬ 
tracted. The values of C for the absorption line spectra as a 
function of aperture size are shown in the left panel of Fig¬ 
ure 9. Also displayed is a guide that indicates the physical 
size probed by a 3 arcsec aperture at various redshifts. 

As Figure 9 shows, the cross-correlation values for the ab¬ 
sorption line spectra (the blue spectra) drop for most galax¬ 
ies as the aperture narrows. This result implies, not surpris¬ 
ingly, that the spectra of the inner regions of some galaxies 
are different than the integrated spectra. The decrease in 
the correlation value at its most extreme is about 50%. In 
Appendix A, we showed that lowering the mean correlation 
value to 59% (S/N = 2.3) reduced the fraction of successful 
classifications to below 70%. For redshifts 2 ; < 0.03 (9000 km 
s^^), the Figure implies that aperture bias may significantly 
distort the classifications of some spiral galaxies. However, 
there are no systematic errors in the classifications of the CL 
sample at low redshifts. This success results partly from the 
domination of the sample by early-type galaxies and partly 
from the classification of a galaxy as late-type regardless of 
its absorption spectrum if it has strong emission lines. 

Several other factors affect the cross-correlation values. For 
example, there are differences in the line profiles between 
the inner and integrated apertures. Line profile differences 
reflect kinematic changes and have a smaller effect on the 



Figure 9: The corrected correlation as a function of aperture 
size for the AB sample (definition of corrected correlation 
is given in Appendix B). The scales within the plots show 
the corresponding physical size of a 3 arcsec aperture for a 
galaxy at the labeled redshift. 


classification than differences in line intensity ratios. For ex¬ 
ample, in cases where rotation causes a bi-modal line profile, 
the classification will be unaffected if both the early-type 
and late-type components are similarly distorted, but the 
cross-correlation value will be affected. Imperfect sky sub¬ 
traction may also contribute to the decrease in the correla¬ 
tion value in at least several cases. Therefore, it is an overes¬ 
timate to attribute all of the decline in the cross-correlation 
values in Figure 8 to differences in the radial distribution of 
stellar populations. 

We treat the emission line spectra (the red spectra) slightly 
differently. We only consider the region of each spectrum 
around the Ha and [N II] lines. The rest of the procedure 
is identical to that used for the absorption lines. The re¬ 
sults are shown in the right panel of Figure 9. Most of the 
galaxy spectra change little between narrow and wide aper¬ 
tures, which indicates that the galaxies in this sample have 
emission lines (or lack them) throughout their disks. 

One notable exception is not well illustrated by the Fig¬ 
ure. In this galaxy, the continuum emission drops below sky 
at a radius at which no line emission is detected, and so 
the galaxy appears to have similar spectral properties at all 
radii. However, the 2-D spectrum reveals strong Ha emis¬ 
sion well beyond the radius where the continuum emission 
is undetectable. This result is a clear example of aperture 
bias. The bias against observing this line emission would 
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occur even if the galaxy had a redshift of several tenths. For 
such a galaxy, the emission lines that would identify it as 
a spiral will not be found in either multifiber or multislit 
data. Although there will always be some fraction of spi¬ 
rals that are misclassified due to this problem, twelve of the 
eighteen spirals we examine have easily detectable emission 
lines, and, of those, only the one has emission exclusively at 
large radius. 

We conclude from Figure 9 that the spectral classification 
of spiral galaxies without detected emission lines and with 
recessional velocities below 15000 km s^^ has about a 50% 
chance of being significantly affected by aperture bias. Com¬ 
bining this percentage with the observation that emission 
lines (affected by aperture bias in one out of 12 galaxies) 


are detected in 70% of the CL sample spirals of morpholog¬ 
ical type Sb or later, we expect that aperture bias may be 
a problem in about 20% of spirals with 2 : < 0.05. Ellipticals 
and lenticulars, which are much more radially homogeneous, 
are even less likely to be affected. Aperture bias does not ap¬ 
pear to be a significant problem in the CL survey, although 
it may explain the 12 morphologically late-type galaxies that 
are spectrally classified as E’s or SO’s. 

We conclude that, for 3 arcsec apertures, aperture bias does 
not seriously affect classification for the vast majority of 
galaxies, even those at redshifts 2 ~ 0.03. These results are 
general trends, and we stress that aperture bias may still be 
very important in certain cases. 



